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Abstract 

The quantum action generated by fermions which are minimally coupled to abelian 
vortex background fields is studied in D = 2 + 1 and D = 3 + 1 Euclidean dimensions. 
We present a detailed analysis of single- and binary-vortex configurations using the 
recently developed method of worldline numerics. The dependence of the fermion- 
induced quantum action on the fermion mass and the magnetic fluxes carried by the 
vortices is studied, and the binary-vortex interaction is computed. Additionally, we 
discuss the chiral condensate generated from a dilute gas of vortices in the intermedi- 
ate fermion mass range for the case D = 3 + 1. As a byproduct, our findings provide 
insight into the validity limits of the derivative expansion, which is the standard 
analytical approach to inhomogeneous backgrounds. 
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1 Introduction 



The dynamics of surfaces plays an important role in many fields of physics ranging from 
solid state physics and chiral quantum field theories |T|, [| to string theory Q]. The 
so-called vortices, d — 2-dimensional surfaces of the ^-dimensional space, are of particular 
interest in view of their phenomenological importance: in solid state physics, freely moving 
vortices give rise to a non-vanishing conductance of high T c -superconductors, thereby lim- 
iting their technical applicability ||. In the context of Yang-Mills theory, it has recently 
been observed in lattice gauge simulations that, in the continuum limit ||, the Yang-Mills 
vacuum is populated by center vortices the core of which can be detected in a certain gauge 
by a projection technique J7|, ||. A random gas of these vortices can grasp the essence of 
quark confinement at zero temperature, and the deconfinement phase transition at finite 



temperatures can be understood as a vortex de-percolation transition ||, |10[. Although 
rigorous estimates signal that the core size of the full (unprojected) field configurations is 
spread over the whole spacetime [ lTJ , an effective vortex model supplemented by a partic- 
ular core size describes certain low-energy observables remarkably well |12|, |13|| . 

In many applications, the surfaces appear as classical background fields (solitons) of 
an underlying field theory. The surface properties as well as the surfaces' interactions are 



determined by thermal and quantum fluctuations, respectively [|14| , [15| . In the context of 
the layered superconductors, the spectrum of the quasi-particles encodes phenomenologi- 
cally relevant information of the vortex background 0, [17]]. In the context of Yang-Mills 
theory, a deep knowledge of the interplay of quarks with vortex-like solitons would help to 
describe hadron properties in the vortex picture. 

The determination of the vortex free energy generically requires the calculation of func- 
tional determinants where the vortex profile enters as a background field. Whereas the 
eigenvalue problem of the single-vortex case with a ^-function profile can be solved ana- 
lytically [[18], a numerical treatment seems inevitable for realistic profiles and multi- vortex 
configurations. One possible technique aims at the numerical integration of a Klein-Gordon 
or Dirac type of equation associated with the differential operator under consideration. For 



example, in |14|, [L5| , this approach has been further developed and the numerical burden 
has been reduced to a quantum mechanical computation of bound-state energies and scat- 
tering phase shifts; in particular, the usually delicate issue of renormalization can advanta- 
geously be inherited from standard perturbation theory in this approach. The disadvantage 
of this approach is that numerical complications increase with the degree of complexity of 
the background field; up to now, only highly symmetric backgrounds have been treated (see 
fH| and |2D[ for the single- vortex case). In this work, we put forward a recently proposed 
numerical technique pi , 22[ which is based on the "string-inspired" worldline method |23) . 
In this formalism, the effective action (functional determinant) is represented in terms of 
first-quantized particle path-integrals, which have turned out to be highly convenient for 
analytical computations involving constant background fields (see |M] for a review). In our 
numerical realization of this formalism, Monte- Carlo techniques can be exploited to esti- 
mate expectation values of background-field dependent operators with respect to scalable 
worldline ensembles, so-called loop clouds. This worldline numerical scheme is formulated 
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in coordinate space, which facilitates a transparent renormalization, since the divergencies 
are associated with local operators in coordinate space and the required counterterms can 
easily be read off. The most important advantage of this formalism is marked by the fact 
that the numerical algorithm can be realized without any reference to the specific choice of 
the background field. In particular, a high degree of symmetry of the background is not re- 
quired at all; on the other hand, if a symmetry is present, it can, of course, be exploited for 
reducing numerical efforts. Most recently, this technique was adapted for worldline loops 



on a cubic lattice, most convenient in the case of lattice gauge simulations [25|. Therein 



it was pointed out that a random walk easily generates loop clouds with the appropriate 
measure required for the worldline approach. 

In this paper, we will explore the free energy of vortex-type solitons, being supported 
by a U(l) gauge field which is minimally coupled to a fermion. We will discuss the case of 
2 + 1 and 3 + 1 spacetime dimensions and restrict ourselves to vortex configurations which 
are static with respect to one or two directions of spacetime, respectively. The free energy 
will be studied for the single and the binary-vortex configuration. In the latter case, we 
will obtain the fermion-induced vortex-vortex interaction. 

The paper is organized as follows: in the next section, the worldline approach to the 
fermion determinant is introduced and the renormalization procedure is discussed in detail 
for the case D = 3 + 1. The vortex configurations under investigation are introduced in 
section [| In section ^, our numerical result for the single- vortex configuration is presented 
and compared with the estimate obtained from existing analytical studies. The binary- 
vortex interaction is presented in section [5[ Fermion condensation due to a dilute gas of 
vortices is addressed in section ^|. Conclusions are left to the final section. 



2 Fermion determinants in the worldline approach 
2.1 Setup 

In this paper, we shall investigate the case of one flavor of a 4-component Dirac fermion in 
D = 2 + 1 and D = 3 + 1 spacetime dimensions. The fermion field is minimally coupled to 
a U(l) gauge field A^x), which is considered as a background field. The fermion-induced 
effective action rf erm , which is a free-energy functional, determines the probabilistic weight 
exp{-r fcrm }. It is given by 

r ferm = - lndet(m - ip) , p := $ + i^, (1) 

where m is the fermion mass, A 1 = A^(x)7 M , and 7^ are the anti-hermitean Euclidean 
7-matrices, 7^ = —7^. As a consequence, pis hermitean. 

Since the Dirac operator exhibits no spectral asymmetry in the 4-component formula- 
tion, the effective action is real, and we obtain in Schwinger propertime regularization 

r fem = ~ r ^e-^Trexpj-^T) , (2) 
1 Jl/A* 1 L J 
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where the scale A acts as a UV regulator. The idea of the worldline approach consists of 
rewriting Eq. (Q) in terms of 1-dimensional path integrals (for a comprehensive review, see 



24]). The worldline representation of Eq. (fj) is given by 



fcrm = 2 (4^ J ^ Q 2 jS^ ^ 4 " l) ; (3) 

where a gauge-field independent constant has been dropped, and we have defined the spino- 
rial counter gauge factor of the Wilson loop including the Pauli term (spin-field coupling), 

W apiQ [A) = ^expjz^ A^x) Xfl drj trP T expQ^ dra^F^. (4) 

Here F^ u (x(t)) is the field strength tensor, and := i[ r y tJ ,, r y u \/2 are hermitean Dirac 
algebra elements. The average (. . . ) x in (|3|) is performed over an ensemble of closed world- 
lines, or loop clouds. A single loop is characterized by its worldline x^(r), r e [0, T] in D 
dimensions. Pt denotes path ordering with respect to the propertime T. The loops are 
centered upon a common average position, 

f T 

Xg := (1/T) / dr x^{t) , ("center of mass") . 

Jo 

The loop ensemble is generated according to the Gaussian weight: 

\ J dr (x\ + . . . + x\ 



cxp 



(5) 



In practice, one greatly reduces the numerical work by generating a loop gas of unit prop- 
ertime T = 1 only, and by producing loop ensembles for a given propertime T ^ 1 by 
rescaling. This numerical method is discussed in some detail in [21], ^ . 



Finally, we point out that the symmetry of the background field can easily be exploited 
to reduce the numerical burden. In our case, the background field is static with respect 
to D — d space-time dimensions, i.e., the gauge field depends only on the coordinates 
xi, . . . , Xd, d < D. Since the weight of the loop clouds is Gaussian, the loop average in D 
dimensions can be evaluated by using d < D dimensional loop clouds only, i.e., 



r[A\) = Vx 1 ... d F[A]e iJo v 1 'J (6) 



where 

f _irT d ( ±2 



In particular, it is sufficient for the examples discussed below to generate a single (large) 
2-dimensional loop cloud to address the properties of static vortices in D = 2 + 1 and 
D — 3 + 1, respectively. 
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2.2 Renormalization 



In the case of D = 2 + 1 dimensions, the propertime integration in Eq. (|3|) is finite at the 
lower bound. This implies that one might safely remove the regulator, i.e., A — > oo, in this 
case. We therefore confine ourselves in this subsection to the case D — 3 + 1, where the 
fermion-induced effective action discussed above is actually the bare action and given by 
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ferm 



2 (Anf 



a xq t=t e 



-m 2 T 



l/A 



2 T 3 



W spin [A] 



(7) 



For small values of the propertime T, the expectation value (W sp i n L4] — l) x can be cal- 
culated analytically for the case of a loop cloud with x being the center of mass. One 
finds 



W spin [A] 
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F ij , u {xq)F ixv (xq) T 2 + ... 



where the ellipsis denotes higher-order invariants that can be formed out of the field 
strength tensor, its dual, and derivatives thereof. Rewriting (|^) as 



ferm 



8tt 2 
1 1 



d 4 x 



OO jrp I 1 

^ e" m2r {W spm [A] - 1 - ^ 2 (* )T 2 



+ 



6 8tt 2 



1/A2 

dT 



-m z T 



l/A 2 T 



J d 4 x F\x ) 



(9) 
(10) 



we observe that the part ((^) of r^ rm is finite if the regulator is removed, A — >• 00. Using 



dT 

1/A 2 T 



-m 2 T 



m 



A 2 



In — - 7E + O — 



77 V 



A 2 J 



(7e is the Euler constant), we neglect irrelevant terms which are suppressed by powers of 
l/A 2 in part © and (0) and obt am 



form 



50 dT _ m 2 T / r I 2, ,- r 2 



48vr 2 



d x Q j — e 
In (^) + 7E 



W spin [A] - 1 - -F 2 (x )T 2 
o 

d 4 x F 2 (x ) . 



(12) 
(13) 



The renormalized effective action r e g is obtained by adding the bare "classical" action 

1 



eff 



ferm ' 



d A x F\x ) 



where g\ is the bare coupling constant, and we enforce the renormalization condition 

1 



1 1 

+ 



9b (A) 



12tt 2 



7e 



(14) 



(15) 
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Here <7R,(/i) denotes the renormalized coupling at a given renormalization point //, and we 
rediscover the QED (3 function /3(g^) = /i<9 M g|(/i) = g^/(6ii 2 ). Inserting fll2"D and (|l5l) into 
(H), the renormalized effective action is given by 



Teff 

Tferm 
r p2 



Tp2 + rf erm 
1 « V 



dT 



-m 2 T 



n T3 
1 



W/ spin L4] - 1 - -F 2 (x )T 
6 



12tt 2 



y d 4 x f 2 (x ) . 



(16) 
(17) 



In these equations, Tp2 denotes the renormalized Maxwell action; here we can impose 
fermion-mass-shell renormalization by choosing /j = m, so that the log term in flT7j ) drops 
out. In addition to the (trivial) classical term (|T?| ), worldline numerics provides us with 
an explicit answer for the renormalized fermion-induced quantum contribution rf erm = 
/ d D xC ierm given by ([16]). 

The renormalization scheme employed here can easily be related to standard renormal- 
ization prescriptions for Feynman amplitudes. The latter are formulated, for instance, by 
considering the renormalized photon propagator D^ v {j>) = P MiV D(p 2 ) in the Landau gauge, 
where P^ u is the transverse projector, and specifying a number R, 



R 



_d_ 

dp 2 



D-\ P 2 



\p 2 =0 j 



which determines the wave function renormalization constant in a particular renormal- 
ization scheme. Since the effective action T e g given above represents nothing but the 
generating functional for one-particle irreducible Green's functions, it is related to the 
renormalized photon propagator by 



D 



-i 



S 2 T eS [A] 



[IV 



SA^x) SA u {y) " 



(19) 



We observe that the part fll6|) does not contribute to the right-hand side of Eq. ([TJ) in the 
limit p 2 — > 0; this leads us to the desired result 



R 



d_ 

dp* 



D-\p 2 



|p 2 =0 



In 



9r(v) 12712 V H 



m 



(20) 



which determines the value of the coupling at any renormalization scale /x for a given value 
olR. 



Vortex interfaces 



The core of a vortex corresponds to a D — 2 dimensional surface of D-dimensional Euclidean 
spacetime. The vortex field which we will discuss below is represented by a U(l) gauge 
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potential A^x). In an idealized case, the gauge potential is singular at the vortex surface 
and can locally be represented by a pure gauge. Considering a closed line C, the vortex 
can be characterized by the holonomy 

exp|i <j> Ap{x) dXfX = exp{i(pL}, (21) 

where is the magnetic flux carried by the surface, and L is the linking number of C with 
the D — 2 dimensional vortex surface. In the case = 7r, the gauge potential describes a 
so-called center vortex in analogy to SU(2) Yang-Mills theory. 

In many physical applications (for instance, in the context of center vortices in Yang- 
Mills theories) one observes that the vortex surface possesses a finite thickness d with 
respect to the directions perpendicular to the flux, implying that the gauge potential 
singularity in the surface is smeared. We point out that the physics of the extended object 
can be quite different from the physics of the idealized vortex. In the latter case, one 
removes the D — 2 dimensional vortex surface from the base manifold. Quantum fields 
are acting in the reduced coordinate space and obey certain constraints at the singular 
subspace. For purposes of illustration consider the case = 2n. One expects that the free 
energy of the vortex is degenerated with the vacuum because the quantum fields of both 
cases are related by a gauge transformation (see [26| for an illustration at the 1-loop level). 



By contrast, in the case of the smeared vortex surface, the quantum fields experience the 
magnetic flux of = 2tt. Therefore one does not expect a degeneracy of the smeared vortex 
configuration with the vacuum. 

In this first investigation, we will concentrate on plane vortex surfaces. For definiteness, 
we consider the gauge potential of a vortex with flux = (fir and core size d of the form 



M X ) = l A2± 2 M 2 fo, -*1> °> > fl=l...D. (22) 



Since the worldline approach is manifestly gauge invariant, any other choice of a gauge- 
equivalent configuration would give the same fermionic interface energy. The corresponding 
field strength is given by 

F12(X) = W + xl + xlf ■ (23) 

In the case of D = 2 + 1 dimensions, the fermion- induced effective action is proportional 
to the extent L t of the space-time in time direction, 

POO 

= EL t = L t 2n / dpp4 3 r ) m (p), (24) 

where p is the radial coordinate of the xy-plane. E can be interpreted as the energy of the 
static vortex soliton. C\H m (p) is a radial energy density, i.e., the amount of energy which 
is stored in a cylindric shell [p, p + dp] . In the case of D = 3 + 1 dimensions and at a given 
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time slice, the vortex core is given by a straight line of length L z . The effective action is 
proportional to the extent L t L z , 

POO 

r fe 4 rm = XU^z = U L z 2tT / dp p C^M , (25) 

Jo 

where x is the string tension of the vortex line in the 3D hypercube. 

4 Quantum energy of the single-vortex configuration 

4.1 The derivative expansion 

The standard analytical approach to effective actions and quantum energies for nonhomo- 
geneous backgrounds is the derivative expansion. Here, the desired answer is expanded in 
terms of a small parameter constructed from derivatives of the background field. In the 
present case, there are two options for the expansion parameter, which can symbolically be 
written as d 2 /m 2 <C 1 or d 2 /B(x) <C 1, i.e., the derivatives of the background field should 
be smaller than the fermion mass or the local field strength. It is remarkable that closed- 
form expressions at the next-to-leading order (NLO) level in D = 2 + 1 and D = 3 + 1 have 
been found in recent years which are applicable in both cases; only one of the conditions 
mentioned above has to be satisfied. The leading-order derivative expansion, which agrees 
with the Euler-Heisenberg effective Lagrangian for constant background fields, is given by 



27, 28 



COO Jrp 

£$L(x) = C dD I —e- m2T [B(x)Tcoth(B(x)T)-l + {c.t.}]. (26) 



The next-to-leading (NLO) correction can be written as [J29], pOlj f] 



.11 /-N _ r, (0iBM) 2 r ^ .--g^ <? 



The sum of both Lagrangians provides us with the NLO derivative expansion of the fermion- 
induced effective Lagrangian. The prefactors Cao,i and exponents a, v and p depend on the 
number D of spacetime dimensions and are summarized in Table [3]. Note that the counter 
term {at.} = — ^(B(x)T) 2 must be added to the integrand in Eq. ( p6|) for D—3+1, as 
discussed in Sect. |2[2| (cf. Eq. (|])); for D — 2 + 1, there is no counterterm, {at.} = 0. 

The expression ( p7|) is valid for any value of the ratio m 2 /B, i.e., either m 2 or B 
has to be large compared with the inverse length scale squared set by the variation of 
the background field. Near the vortex core, this length scale is naturally given by the 
core size d. For a flux (p = 0(1) considered in the following, we find from Eq. (p3|) that 
B(x)d 2 = 0(1) near the core. Simultaneously assuming small masses, md < 1, neither 

1 A representation of the integral in terms of the Hurwitz Zeta function as well as asymptotic expansions 
have been found; see, for example, p0|. 
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D=2+l 



£>=3+l 



Cdo 


i 

47I-3/ 2 


l 

8vr 2 


Cdi 


4(47t) 3 / 2 


1 

(8tt)2 


a 


1.5 


1 


V 


2.5 


3 


P 


0.5 


1 



Table 1: Prefactors and exponents of the derivative expansion. 



of the possible derivative expansion parameters is small. The NLO Lagrangian reaches 
its validity limit, and one cannot expect reliable results from the derivative expansion. 
Below, it will turn out that the quality of the NLO approximation strongly depends on 
the number of spacetime dimensions. By contrast, for large masses md > 1 or for large 
radial distances from the vortex core p > d , we expect the NLO Lagrangian to give a 
reasonable answer, since one of the possible expansion parameters is small. Below, we will 
use the regime of large masses (or large p) to gain insight into the range of parameters to 
which worldline numerics is reliably applicable (see [21j and [^5| for a first comparison of 
worldline numerics to the constant-field case). 



4.2 Worldline numerics: D = 2 + 1 



Closely following the procedure outlined in |21j, we generate a loop ensemble with n e 
loops, where each loop is represented by a set of N spacetime points with coordinates 
x\. . .%n- The loop ensembles are generated with appropriate measure (see discussion in 
subsection PTT| ). We have used 10.000 and 20.000 dummy sweeps and found that for a 
proper thermalization, 10.000 dummy sweeps are sufficient for the accuracy achieved in 
the results shown below. Furthermore, we have used N e 100, 150, 200 points specifying 
each loop. Whereas iV = 100 seems sufficient for our purposes in the (2+l)-dimensional 
case, we take iV = 200 to generate high-precision data for the case D = 3 + 1. Generically, 
we average over n e = 1000 loop ensembles in the (2+l)-dimensional case, and exceptionally 
use n e = 16000 for certain applications in D = 3 + 1. Any dimensionful quantity quoted 
in the following is calculated in the simulation in units of vortex thickness d. 

In order to avoid a violation of gauge invariance in the numerical computation, we have 
to deal with a subtlety concerning the discretization of the integrals along the worldlines: 
a particular loop in our ensemble is considered to be a polygon with straight lines C« 
connecting the points Xi and Xi + \. In order to evaluate the holonomy in Eq. (gj), we 
consider the "infinitesimal" part 



Ci 



expii / Afj,(x) x^ dr > (28) 



9 



2 2 

m =l/d 
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Figure 1: Effective Lagrangian as a function of the radial distance p to the vortex core for 
the case m—1, tp—1, D—2+1. The NLO derivative expansion (solid line) is compared with 
the numerical computation (squares with error bars). 

for each Ci separately and evaluate the integral analytically using the vortex profile under 
consideration. This procedure guarantees that our numerical result is invariant under 
gauge transformations of the background gauge field A^x) for any number N of points 
defining the polygons. Moreover, the flux enclosed by the polygon is exactly taken into 
account as desired (and not only within discretization errors). Of course, this procedure 
can be generalized to arbitrary vortex profiles for which the "infinitesimal" integrations 
along the Cj's can be performed numerically; thereby, the properties mentioned above can 
be preserved to any numerically desired precision. However, we should stress that the use 
of a smooth gauge (e.g., covariant gauges) for the background field is recommended in this 
case; this facilitates a fast convergence of the numerical integration. 

In order to test our numerical approach, we calculate the effective Lagrangian Cf erra (p) 
for large values of p, p ^> d, where the derivative expansion is expected to give reliable 
results. The result of the numerical worldline approach is compared with the derivative 
expansion in Fig. [I] for D = 2 + 1. The agreement between the two curves is satisfactory; 
in particular, the numerical approach is able to compute £f er m over a range of many orders 
of magnitude. In the region close to the core of the vortex, i.e., p ~ d, the gradients of 
the background field are as large as the field itself, so that the reliability of the derivative 
expansion now depends on the value of the mass. Our numerical findings for this regime 
are shown in Fig. |2| for m? = 1 (left panel) and m 2 = 5 (right panel). For larger masses, we 
observe a good qualitative and a reasonable quantitative agreement. But even for the small 
mass value m 2 = 1, there is at least qualitative agreement between the numerical result and 
the derivative expansion, indicating that the applicability of the derivative expansion can 
be pushed to its formal validity limit in D = 2 + 1. This observation is also supported by 
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m=l/d 



2 2 

m =5/d~ 




the fact that the NLO term ( p7|) is only a small correction to the zeroth-order result (|26|). 
Moreover, we expect that the (up to now unknown) NNLO correction, which is sensitive to 
the curvature of the field strength, improves the result near the vortex core. In this sense, it 
is reassuring to observe that the numerical result agrees with the NLO derivative expansion 
precisely at the turning point of the curve at p ~ 0.5, because the NNLO correction must 
vanish here. We should finally stress that for even smaller masses m < 1, the discrepancy 
between our numerical result and the derivative expansion increases, so that the derivative 
expansion should be abandoned here. 

Let us now examine the fermion- induced quantum energy E of the vortex soliton as 
defined in Eq. (H, 

/•OO 

E = r®JL t = 2n dp P C lcim (p), (29) 
Jo 

which we obtain by numerically integrating the effective Lagrangian. Our result for this 
energy is shown in Fig. |3] as function of the fermion mass m in units of the vortex thickness 
d. Since fermion fluctuations are suppressed with increasing mass, the quantum energy 
decreases with increasing m, and vanishes in the large mass limit. For phenomenological 
purposes, it is important to notice that the quantum energy is positive. This implies that 
potential effective models for vortex dynamics have to account for the fact that vortex 
nucleation is suppressed by the fermion-induced effective action in D = 2 + 1 dimensions. 

Furthermore, let us consider the variation of the quantum energy with respect to the flux 
4>, which is carried by the vortex. Our numerical result is shown in Fig. (right panel). We 
find that the energy is monotonically increasing with the flux 0. For <fi = <£>7r=27r, the vortex 
configuration approaches asymptotically (p — > oo) a pure gauge. As in the instanton case, 
the energy is nonvanishing due to the finite extension of the vortex core. It is interesting to 
compare our result for E((p) with a general result for fermionic determinants described by 
M. Fry in |HJ]; therein a lower bound has been derived for unidirectional magnetic fields in 
D — 2 + 1, which translates into an upper bound E\, for the quantum energy of our vortex 
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0.6 




Figure 3: Quantum energy as a function of the mass m 2 for (p=l (left panel) and as a 
function of the flux (p carried by the vortex for m = 1 (right panel) in comparison with 
Fry's upper bound given in Eq. (p0|). 



configuration given by 

E(<p) < E h {<p), (30) 
E h (cp) = Q j i 2 - 3cp - 2^/1 + cp + + cp + 3^ Arsinh(V^) 

for the case md = 1 and Dirac 4-component spinors. For other values of md, this formula 
receives a total factor of (md) 3 and the flux has to be replaced by (p — > ip/(md) 2 . As shown 
in Fig. |3|, our numerical result lies well within this bound. More remarkable is the fact that 
the functional dependence of our result agrees with the bound within the error bars, if the 
bound is scaled by a factor of roughly 0.65. 

As a further check, we have compared all our above-mentioned results with those of 
H| , where the single- vortex case with profile functions different from ours was considered 
within the phase-shift approach. We find good agreement within the error bars except for 
a global factor of 2 by which the result of [B5| for E(ip) is larger. 



Let us finally point to a problematic feature of our present approach. Note that our 
statistical error bars increase in the regime of small fermion masses in Fig. | (left panel), 
which has the following origin: the Pauli term in Eq. (HI), 



trP T exp ( - / dra^F^ 



o 



can favor large loops (T ^> d 2 ) for the present unidirectional magnetic field. On the other 
hand, the holonomy factor, 

exp iil (x) dr 



o 



changes its sign rapidly in this case when the loop under consideration is slightly deformed. 
This implies that a finite result in the small-mass regime arises after subtle cancellations, 
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Figure 4: Effective Lagrangian for m 2 = 1 (left panel) and m 2 = 3 (right panel), (p—1, 
D = 3 + l. 



and is therefore hardly accessible to the present numerical formulation. For finite val- 
ues of the mass, potentially large Pauli term contributions are suppressed by the factor 
exp(— m 2 T), which solves the numerical problem for large loops (large T).0 However, this 
cancellation problem can become serious for the approach to the chiral limit; here, a solu- 
tion to the problem has to be implemented in the algorithm on the analytical level rather 
than by brute-force numerical means. 

4.3 Worldline numerics: D = 3 + 1 

Similarly to the previous studies of the D = 2 + 1 case, we investigate the effective La- 
grangian £ferm as a function of p and compare the result with the one obtained from the 
derivative expansion (|26|j27|) (see Fig. f|). Contrary to D = 2 + 1, we observe that unless 
p ^> d, the leading order (LO) of the derivative expansion falls far too short of reproducing 
the numerical result for m 2 = 0(1). For instance, for m 2 = 1, we find an order-of- 
magnitude difference at p = (Fig. |] (left panel)). Moreover, the NLO contribution of the 
derivative expansion also exceeds the leading order result by almost an order of magnitude 
for p « d. This signals the break-down of the derivative expansion for the case of the field, 
its gradient and the mass being all of the same order. Contrary to the case of D — 2 + 1, 
the applicability of the NLO derivative expansion cannot be pushed to its formal validity 
limits md ~ 1. Even for larger values of the mass, m 2 = 3 (right panel), there is only a 
minor improvement of the quality of the NLO derivative expansion. Of course, the NNLO 
contribution could, in principle, improve the results of the derivative expansion, but this 
would only emphasize the fact that there is no clear hierarchy from term to term in the 

2 In the constant-field case, this cancellation problem already occurs for m 2 < B as observed in p^ |. In 
the present vortex case, we can afford much smaller masses. In general, the smaller the coherence length of 
the field, the less serious the cancellation problem. The usually considered homogeneous and unidirectional 
field configurations are therefore rather pathological in this respect. 
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Figure 5: Effective action as a function of the dimensionless quantity m 2 d 2 , <p=l, D=3+l 
(left panel), and as function of the flux, m 2 d 2 = 1 (right panel). 

derivative expansion. 

We believe that the striking difference to the D = 2 + 1 dimensional case is indeed 
remarkable and points to a deeper reason in terms of a renormalization effect. To illus- 
trate this, we note that the quantity (Ws P i n — 1) occurring in the propertime integrand 
is positive for the vortex background, whereas the counterterm — ^B 2 (x)T 2 is negative. 
Since the effective Lagrangian is largely negative as seen in Fig. ^, it is mainly driven by 
the counterterm. Now the leading-order derivative expansion obviously overestimates the 
value of (W sp i n — 1) near the vortex core, since it is a local expansion. The true value as 
seen in the numerical computation is much smaller because it is a nonlocal average over the 
extended loop cloud that also "feels" the much weaker field at a radial distance from the 
core. The final value of the total effective Lagrangian at a point x therefore results from a 
nontrivial interplay between nonlocal (and nonlinear) vacuum polarization (~ (Wspin ~~ 1)) 
and a local definition of the coupling giving rise to a local counterterm. In regions where 
the background field varies rapidly, such as the near vortex core in our case, this interplay 
can lead to an order-of-magnitude enhancement of the effective Lagrangian as compared 
with the constant-field approximation (leading-order derivative expansion). In our opinion, 
this phenomenom clearly deserves further investigation. 

Returning to our numerical study of the one- vortex background, we calculate the string 
tension x as defined in Eq. (^) as a function of the fermion mass m and plot it in Fig. |5| 
for ip = 1. The negative values of x show that the fermion-induced effective action Ff erm 
favors the nucleation of vortices. Since the modulus of this effective action increases if 
the vortex thickness d is decreased, the fermionic part of the vortex action supports the 
existence of thin vortices. These results are in contrast to those of the case D = 2 + 1, 
where the effective action turned out to be positive (see Fig. [$]). This sign difference is 
again related to renormalization effects mediated by the counterterm. 
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anti-aligned flux 

E = 0.9 2" E , n=0..9 




aligned flux 

E = 0.9 2" E , n=0..9 




Figure 6: Effective Lagrangian L\J m (x,y) as function of the xy-plane for the two vortex 
configuration: parallel (right panel) and anti-parallel (left panel) orientation of the flux, 
m 2 d 2 = 0.5, <p=l, D=2+l. 



5 Binary-vortex interactions 

Since we are investigating the case of Abelian gauge configurations, the binary-vortex 
configuration is given by the superposition of two single vortex gauge fields A^(x) (p2|), 



Af{x) 



A, 



r l - ) + A 



l_ 

m x + 2 



(31) 



where I denotes the vortex distance. Below, we will study the case of the so-called center 
vortices the flux of which is given by <p = 1. The relative sign between the gauge fields on 
the right-hand side of Eq. ( |3"T| ) corresponds to the relative orientation of the fluxes: the 
plus sign corresponds to an equal orientation of the flux in each vortex, while the minus 
sign signals an opposite orientation. 

Figure | shows the lines of equal effective Lagrangian (x, y) in the xy-plane which 
is perpendicular to the vortex fluxes. The vortices are located at the x axis at a distance 
/ = 3d. It is straightforward (but computer time consuming) to integrate the effective 



action over the iq/-plane in order to derive the quantum energy E (|24]) of the binary-vortex 
configuration. The result is shown in Fig. |7| for the case D = 2 + 1. For large distances 
/ ^> d, the quantum energy approaches twice the value of a single vortex. For / = 0, the 
vortices fall on top of each other. If the fluxes of the vortices are oppositely oriented, the 
vortices annihilate each other, and the quantum energy of the configuration vanishes. If 
the vortex fluxes are equally oriented, the configuration is equivalent to the single vortex 
configuration with flux if = 2. Since the quantum energy is roughly proportional to if 2 
(see Fig. Q), the vortex configuration with flux (p = 2 possesses a higher energy than twice 
the energy of a single vortex, carrying flux if = 1. Hence, vortices with an equal flux 
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Figure 7: The interaction of two parallel vortex lines in D = 2 + 1, ip = 1, m cr = 1. 



orientation repel each other in D = 2 + 1, while vortices with oppositely oriented flux 
attract each other. The same line of argument applies to the case D — 3 + 1. Since the 
fermionic contribution to the effective action is negative in this case, the fermion- induced 
force is attractive (repulsive) for equally (oppositely) oriented vortices, contrary to the 
D = 2 + 1 case. 



6 Fermion condensation in a vortex gas 

The fermion condensate can be related to the effective action T c g (JTTJ) by differentiation 
with respect to the fermion mass m, 

[ A (#) = (32) 
J am 

Let us concentrate on D = 3 + 1 dimensions in the following. Here the condensate reads 
in worldline representation: 

m J d*x = £^ / d 4 x Q j^e- m2T (w spin [A] ~ l) , (33) 

where we have first performed the mass differentiation at an arbitrary renormalization 
point and then implicitly chosen on-shell renormalization. 

In the present section, we intend to estimate the quark condensate which is generated 
by a dilute gas of vortices of flux ip = 1 in D = 3+1. For this purpose, we first calculate the 
contribution of the single- vortex configuration with the help of Eq. (|33|) . Due to translation 
invariance of our single-vortex background, the string tension x (p5|) is provided in units 



16 



of the vortex core size only. This allows us to define the dimensionless function cq by 



m J d 4 x = L t L z ^ c (m 2 d 2 ^ , (34) 

where Cq can directly be obtained from worldline numerics, 

c,{ m 2 d 2 ) = 2, f^dQ t |°°|e-W 2 ^W 8pin [^]-l) • (35) 

Here all dimensionful quantities are scaled in units of the core size d, e.g., T = T/d 2 , and 
p = \/ x 2 + y 2 measures the radial distance from the vortex core. In the following, we 
are interested in a dilute gas of vortices which are static and aligned in the z direction, 
but intersect the xy plane at random locations with random fluxes if = ±1. In this 
plane the vortex gas can be characterized by a planar vortex area density py, which is the 
average number of vortices per xy unit area. The total number of vortices within the 4- 
dimensional spacetime is given by Ny = py L x L y . We expect the dilute-gas approximation 
to give reasonable results, if the average distance between two neighboring vortices is at 
least > 2d, where the fermion-induced vortex interactions become small (see, e.g., Fig. |7| for 
the D = 2 + 1 dimensional analogue). In other words, there should be less than one vortex 
per core-size area, py (nd 2 ) < 1. In this dilute-gas approximation, the fermion condensate 
averaged over spacetime volume V is therefore given by 

« V Nv J = 4^2 m Pv c °( m2 d2 ) ■ ( 36 ) 

In the limit of large fermion masses (m ^> 1/d), Eq. (|55D can be studied analytically with 
the aid of the heat-kernel expansion (H), and we find a contribution to the condensate 
which is subleading in 1/m: 

(H) *2^^J F 2 (x ) + o(l/m 2 ^ . (37) 



Here, we recover the familiar result |32| that in the large-m limit the fermion condensate is 
proportional to the field strength squared of the background field (i.e., "gluon condensate" 
in a QCD-like language). 

We investigate numerically the interesting regime of small masses m, where the heat- 
kernel expansion breaks down. The representation (|35| ) is highly convenient for this pur- 



pose. However, as already pointed out, the loop average (...) in Eq. (|34D is plagued from 



a severe cancellation problem which in the present case of the loop parameters employed 
here limits the region of validity to T/d 2 < 100. In order to ensure this limit, we confine 
ourselves to the mass range m 2 T max > 10, i.e., m 2 > 0.1/d 2 . The investigation of the 
small mass regime (m 2 < 0.1/d 2 ) clearly needs further study. Work in this direction is in 
progress . 
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Figure 8: Co (0) (left panel) and the fermion condensate (right panel) as function of the 
fermion mass m, (p—1, D=3+l. 



Our numerical result for cq is presented in Fig. For small masses, i.e., 0.1 < m 2 d 2 < 
0.5, we find that co(m 2 <i 2 ) ~ l/(md), which indicates that the condensate approaches a 
plateau in this regime according to Eq. (|36|), 



m 2 d 2 =O(0.1) 



4?T 



m p Co I m 



,2 j2 



d< 



m 2 d 2 =O(0.1) 



0.058 



1 



47T 2 



pi d 



(38) 



as depicted in Fig. 
the chiral limit m 



right panel). However, we do not believe that this result extends to 
0. Guided by the case of a constant background field in D = 3 + 1 
35| , we expect that the chiral condensate vanishes according to (ipip) — > mlnm for m — ► 0. 
As seen from Fig. |^, our numerical result for Co does not discriminate between the lnm 
and the desired 1/m behavior in the small-mass regime. 

Let us stress the two main findings of this subsection: first, the vortex-gas- induced 
condensate is characterized by a comparably low value of Co for all mass values depicted 
above and an onset of the plateau value for comparably large masses m 2 d 2 ~ 0.1; we 
point out that it is not possible to fit the form of the condensate depicted above within 
a derivative expansion even qualitatively by arbitrarily varying B. Second, in the context 
of a center vortex model for low-energy QCD, the current quark masses are in fact finite, 
though small. The present result indicates that contributions to the chiral condensate have 
to be expected from an interaction of the massive quarks with the gauge vacuum (modeled 
by a random vortex background). As a first estimate, we may insert parameter values 
known from lattice calculations. In the case of an SU(2) gauge theory, the density of the 
center vortices is roughly given by pv = 3.6/fm 2 . The planar vortex correlation function 
was also studied in SU(2) lattice gauge theory fl34"f . One finds an exponential decrease of the 
correlation function. This allows for a definition of a vortex thickness of d ~ 0.3 . . . 0.4 fm. 
Since pynd 2 = O(0.1) in this case, the dilute-gas approximation should be applicable. 
From Eq. (03), we find for the vortex-induced condensate (qq) ~ [50MeV] 3 , for masses 
m 2 d 2 = 0.1, i.e., m ~ 100 MeV. As discussed above, the present algorithm cannot treat 
the case of light mass values, 



m 



5 ... 10 MeV. 
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The vortices considered here obviously possess a trivial topology, i.e., they do not 
contribute to the topological charge. The origin of a chiral condensate here is entirely 
relegated to the field strength carried by the vortex cores. In a fully non-Abelian QCD-like 
context, topologically non-trivial vortex configurations (due to the presence of low-lying 
fermionic modes of the Dirac operator) as well as multi-gluon exchange interactions can 
be expected to provide for a drastic enhancement of the condensate. 



7 Conclusions 

The fermion-induced quantum action of Abelian vortex configurations has been studied in 



the case of D = 2 + 1 and D = 3 + 1 dimensions using worldline numerics [^T], [22], p5| . 
The quantum action of a single-vortex configuration is characterized by the fermion mass 
m, the vortex thickness d and the flux (p carried by the vortex. Our numerical approach 
has been successfully tested in the parameter regime where the derivative expansion is 
expected to provide reliable results: in the large-mass regime md > 1 or for strong-field 
suppression of the inhomogeneities d 2 /B <C 1. 

Compared with our numerical results in D = 2 + 1, the derivative expansion provides 
for a reasonable approximation to the quantum energy of a single vortex configuration even 
for smaller masses, md < 1. The next-to-leading order only adds a small correction to the 
leading-order result. By contrast, in D = 3 + 1, the derivative expansion is insufficient 
even for comparably large masses, md = 0(5). Only for very large masses, md ^> 1, can 
the derivative expansion be trusted. We have argued that this disparity between D = 2 + 1 
and D = 3 + 1 arises from renormalization effects in the latter case, which can occur near 
regions where the background varies rapidly. In particular, these renormalization effects 
can lead to an effective enhancement of the quantum action. Further investigations are 
planned to settle this issue. 

From a physical point of view, we find that the quantum action is positive in D = 2 + 1, 
implying that the presence of vortices is suppressed, and large vortex core sizes d are 
preferred. In D = 3 + 1, the properly renormalized quantum action turns out to be 
negative: the nucleation of thin (d — > 0) vortices in D = 3 + 1 is supported by the fermion 
induced quantum action. 

Subsequently, the binary- vortex interaction was investigated. We found that the fermion- 
induced interaction favors vortices with an opposite orientation of the fluxes in D — 2 + 1, 
while in D — 3 + 1, a unique orientation of the fluxes is preferred. 

It should be stressed that these statements refer to and are derived from the fermion- 
induced action. In a pure QED context, the classical action has to be taken into account. 
The latter will dominate the fermion-induced action by far at weak coupling, reflecting 
the usual hierarchy between classical and quantum-induced nonlinear electrodynamics. 
Formally, the classical and quantum action can become comparable in magnitude for ex- 
ponentially strong fields in the vortex core (exponentially small d for fixed flux); however, 
this is nothing but a manifestation of the Landau pole of QED and thus should be rated 
as unphysical. 
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Finally, the fermion condensate which is generated by a dilute gas of vortices was 
studied as function of the fermion mass md > 0.1 in D = 3 + 1. We found that the 
condensate decreases like 1/m for large values of the mass and that it reaches a plateau 
for md < 0.5. Similarly to the case of a constant magnetic field, we expect that the 
condensate vanishes like m lnmin the chiral limit m — > 0. However, our results point to an 
interesting phenomenon: although chiral symmetry breaking might be tied to topological 
properties of the background fields not covered by the present considerations, the field 
strength carried by the vortex cores enhances the chiral condensate in the intermediate 
fermion mass regime. In a center vortex model of low-energy QCD, this effect may lead 
to a nonnegligible contribution to the condensate. We regret that our worldline numerical 
algorithm in its present form cannot address the small-mass regime md < 0.1 due to severe 
cancellations. The study of the chiral limit in general and chiral symmetry breaking by 
vortex background fields in particular requires improved algorithms and is left to future 
work 
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